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Abstract 

We investigate the relativistic model of superconductivity in (2+l)-dimension. We employ the 
massless Gross-Neveu model at finite temperature and density, to study the superconductivity 
and superconducting instability. Our investigation is related to the superconductivity in (2+1)- 
dimensional two-band systems like MgE>2 or intercalated graphite. 
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I. INTRODUCTION 



Condensed matter physics in low-dimension is an interesting research area both experi- 
mentally and theoretically [1]. For example, the discovery of the high-T c superconductors 
made a large impact on condensed matter physics, and in such systems, the superconductiv- 
ity is considered to occur in a quasi- two-space-dimensional copper oxide Cu0 2 plane [2-7]. 
The recent discovery of the superconductivity in MgB 2 was also a quite important event [8] , 
and people expect that a plane constructed by B atoms plays the main role in the super- 
conductivity. Graphite and carbon nanotube also attract our attention, and they are plane 
systems [9] . Therefore, the superconductivity in two-dimensional systems is of prime impor- 
tance, and well recognized in condensed matter physics. In the theoretical side, Aoki et al. 
studied the effect of the dimensionality [10], the effect of band structures [11], the multiband 
effect ( two-band model [12], four-band model [13] ) and the effect of shapes of the Fermi 
surfaces [14], in superconductivity of several systems ( see also Ref. 15 for two-band model 
and Ref. 16 for three-band model ). The keywords of recent theoretical investigation into 
superconductivity are "two-space-dimensional" and "two-band." 

The two-band models of superconductivity [17-28] have the origin in the papers of Suhl, 
Matthias and Walker [17], and Kondo [18,19], and they were applied to several systems 
under various situations. Both MgB 2 and graphite have honeycomb lattice structure, and 
essentially they are two-band systems. Quite recently, a theoretical study of two-band 
superconductivity in MgB 2 appeared [22-24]. The experimental evidence for two-band su- 
perconductivity in MgB 2 was also obtained [29]. Some relations between superconductiv- 
ity and excitonic state ( exciton condensed state ) in two-band models were also exam- 
ined [25,26]. Low-energy effective theories for two-band superconductivity similar to the 
Ginzburg-Landau model were proposed [27,28]. Some theoreticians consider that the cop- 
per oxide high-T c superconductor can also be described by a two-band theory ( the d-p 
model ) [3,4]. Today, much attentions are paid for two-band superconductivity. There are 
various two-band models for superconductivity. For example, Suhl et al. used the model 
which has only an attractive interaction between particles in a two-band system [17]. Yamaji 
discussed a pairing problem by an interband polarization function arised from a repulsive 
interaction [20]. Kondo found that a kind of two-band effect enhances the superconductiv- 
ity [18]. However, until now, it is not clear how much is the contribution of the lower band. 
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There is no quantitative understanding about the strength of the two-band effect, especially 
the lower-band effect. 

Relativistic fermion often appears in condensed matter systems [30-36] . Semenoff studied 
the (2+l)-dimensional ( 2 for space and 1 for time ) relativistic fermion model in graphite. 
Based on the character of (2+l)-dimensional relativistic quantum field theory [37-41], he 
discussed the anomaly ( the Chern-Simons term and the fractional fermion number ) [30]. 
As a consequence of the honeycomb lattice structure of graphite, its band structure has 
two degeneracy points in the first Brillouin zone ( two conical intersections between the 
upper band and lower band ). Then the relativistic fermion model is obtained in the linear 
dispersion approximation [30-34]: 

e = VFl °( 1 1 Px + 1 2 Py ) (1) 

( as illustrated in Figs. 1 and 2 ). Here, vf is the Fermi velocity, &ndp x andp y are momentum 
operators. By using the Dirac fermion model with a Coulomb repulsion, Gonzalez et al. 
concluded that the superconductivity can emerge in graphite [32]. Shankar derived a Dirac 
fermion model in (2+l)-dimensional doped antiferromagnets [35]. Because those relativistic 
models have two bands ( positive and negative energy states ), they can provide powerful 
techniques to study the low-energy and long-wavelength properties of (2+1) -dimensional 
two-band systems. 

With these considerations given above, we investigate the relativistic model of two-band 
superconductivity in (2+l)-dimension. The purposes of this paper is to examine a two- 
band effect in superconductivity in the (2+l)-dimensional system. We concentrate on the 
examination of the effect of the lower band, by extracting the strength of its contribution in 
the superconductivity. In real substances, our theory can be applied to systems which have 
two- space-dimensional honeycomb lattice structure like MgB 2 , electron-doped-graphite and 
graphite intercalation compounds ( LiC 6 , KC 8 , etc. ), or to systems of the kagome lattice 
structure [42]. Because real substances are not exactly (2+l)-dimensional, we do not consider 
the Kosterlitz-Thouless transition in two-space-dimensional superconductivity seriously [43]. 

This paper is organized as follows. In Sec. II, we introduce the Gross- Neveu model as our 
model Lagrangian, and discuss its characteristic aspects. By using this model Lagrangian, 
we study the superconductivity, namely, one of the phenomena of dynamical ?7(l)-gauge- 
symmetry breaking in (2+l)-dimension. In Sec. Ill, the Gor'kov formalism [44] for a 
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contact attractive interaction in our theory is presented. In Sec. IV, the group-theoretical 
consideration of the mean fields ( gap functions ) is provided. In Sec. V, by using the Gor'kov 
formalism, the gap equations in our theory are derived, and they are solved numerically. In 
1965, Kohn and Luttinger proved that, whether a two-body interaction is attractive or 
repulsive, there is a Cooper instability in an interacting many-fermion system ( the Kohn- 
Luttinger effect ) [45]. Shankar and Chubukov independently proved the existence of the 
Kohn-Luttinger ( KL ) effect in a two-space-dimensional system [46,47]. In Sec. VI, for 
the pairing problem in the case of a repulsive interaction, the KL effect in our theory is 
examined by using the Bethe-Salpeter ( BS ) formalism. Finally in Sec VII, we give the 
conclusion of this work, with further possible investigation in our theory. 



II. THE MODEL LAGRANGIAN AND HAMILTON! AN 

In this section, we discuss the character of the model Lagrangian to study the (2+1)- 
dimensional relativistic superconductivity. We take the following Lagrangian for the starting 
point: 

C(x) = ^{x)i^d^{x) + ^{ X )^ x )f. (2) 

This is the (2+l)-dimensional Gross- Neveu model [48-52]. The first term is the kinetic term 
of the Dirac field, ip and ip are the two-component relativistic spinors describing the Dirac 
fields. Here, we do not give a mass to the fermion. It is well known in (2+l)-dimensional 
relativistic field theory [37-41], the Dirac mass term m^ty violates both the parity and the 
time reversal symmetries. As discussed in the introduction, we treat the system which has a 
degeneracy point at zero momentum. We have to consider the massless case. In Eq. (2), we 
introduce the four-body contact interaction at the same spacetime point. Here, we consider 
one of the simplest relativistic interaction which may generate a superconductivity. When 
G > 0, it will give an attractive interaction. We set aside the question of the origin of the 
attractive interaction, and regard the Lagrangian (2) as a phenomenological one. On the 
other hand, when G < 0, it will become a similar interaction to the on-site repulsion of 
the Hubbard model in its continuum limit. If we employ naive power counting analysis, we 
find that the theory is unrenormalizable and we have to introduce a cutoff. Until now, the 
(2+l)-dimensional Gross-Neveu model is studied extensively, especially in the context of the 
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dynamical chiral symmetry breaking [48-52]. Rosenstein et al. treated the model by using 
the 1/N expansion ( here, N is a number of flavor ) [49]. The method of 1/N expansion, 
a kind of mean-field theory, treats the sum of an infinite number of Feynman diagrams 
similar to the Hartree-Fock theory. Rosenstein et al. also gave a proof of renormalizability 
of the model under the 1/N expansion. On the other hand, we study a subset of diagrams 
different from that of Rosenstein et al. In fact, Rosenstein et al. investigated the chiral 
symmetry breaking arised from fermion-antifermion condensate (ipip) 7^ 0, while we treat a 
fermion-fermion pairing problem, (ipip) ^ 0. In other words, Rosenstein et al. studied the 
dynamical generation of Dirac mass, while we study the dynamical generation of Majorana- 
type mass. In this paper, we will not study the possibility of renormalizability of our 
theory, and treat ultraviolet divergences with a simple cutoff scheme. The cutoff indicates 
the momentum range appropriate to the linear dispersion approximation. Our Lagrangian 
should be regarded as a model in this range of momentum space. For the purpose of this 
paper, we do not have to introduce the constant for the Fermi velocity Vf in the model. 
Equation (2) itself has symmetries of Poincare invariance, L r (l)-gauge invariance, charge 
conjugation invariance, spatial inversion and time reversal invariance. 
For the Clifford algebra [37,38] of 7-matrices, using the Pauli matrices: 

7° = ^, ^=za\ 7 2 = i<7 2 , (3) 

( here the chirality rj = |tr7°7 1 7 2 = +1 ) we obtain the next relations: 

h", y} = 2«r , 7 M = 7V7 , w = <r - ^ V7 p . (4) 

Here, e 012 = 1. We take the metric convention as g^ v = diag(l, — 1, — 1). The charge 
conjugation matrix is given as 

c~Vc = -Y T , C^C = 1. (5) 

For the Poincare algebra, the generator of the Lorentz transformation satisfies the 5*0(2, 1) 
algebra: 

[f,f] = -ie^ x g Xp f, (6) 

where j° is the generator of two-dimensional rotation as the U(l) phase transformation, 
while j 1 and j 2 are the boost operators. Especially the representation on the Dirac field is 



given as 



(7) 



then ip transforms as 



(8) 



The Hamiltonian of our problem becomes 



H 



d 2 xi)(-tf ■ V - tV)^ - 77 




(9) 



Here, we introduce the chemical potential [i. In the relativistic field theory, it descibes the 
finite density at fi ^ as the conjugate of the particle number minus antiparticle number. 
When we consider a crystal by the model, \i becomes the conjugate of the electron number in 
upper band minus hole number in lower band. \i determines the position of the Fermi level 
measured from the degeneracy point. The system has a Fermi circle at ji ^ 0. Throughout 
this study, we completely neglect the temperature dependence of /x, and we treat \x as a 
parameter introduced from the outside of the system. Thus we set /i = ep ( Fermi energy 
). The Fermi level of graphite locates on the degeneracy point. Hence, in the case of 
intercalated graphite, \i determines the electron doping concentration. 

Before closing this section, we would like to comment on the reason why we use the 
theory of the two-component Dirac field. In the theory of (2+l)-dimension, we have a 
choice between a two-component Dirac field and a four-component Dirac field. In fact, 
the dispersion of the honeycomb lattice has two degeneracy points. Gonzalez used a four- 
component theory to take into account the effect of the interaction between the two Fermi 
points in graphite [32]. However, the purpose of this paper is to extract the two-band 
effect, especially the lower-band effect. For this purpose, we have to construct the model as 
simple as possible. The criteria of our model are given as follows: (1) It will give a conical 
( relativistic ) dispersion, (2) it has an interaction which may arise a superconductivity. 
From these criteria, we consider the Lagrangian (2) as the simplest model. If we use a four- 
component theory of superconductivity, it gives a problem of the effect of the interaction 
between two distinct Fermi circle. We regard that the effect is essentially different from the 
lower band effect which is studied in the context of this paper. Hence we use a two-component 
theory, and treat the problem of a single conical-dispersion system. The operator 7 M p M has 
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two linearly independent solutions; one for positive energy state and another for negative 
energy state. Therefore, though ip is a two-component spinor, there is no degeneracy of a 
spin-degree of freedom. Both the upper and lower bands have no spin degeneracy. We treat 
a kind of spinless model, and it is enough for our purpose of this paper. 



III. GOR'KOV FORMALISM 

In this section, we derive the Gor'kov theory [44] for pairing problem under the attractive 
interaction G > in the Lagrangian (2). The formalism given in this section is parallel with 
the (3+l)-dimensional theory [53-55]. The field equations are obtained from the Lagrangian 
(2) by the action principle: 

i^d^ + Gf^)^ (10) 

-id^ + G^)^. (11) 

For the derivation of the Gor'kov formalism, we introduce various propagators. We use the 
4-component Nambu notation [56]: 

( t ( rl ) ' = (^)^))' 

\ V (x) J 

where T means the transposition. The definition of the one-particle propagator is 

G(x,y) = -i{T*(x)*(y)) 








dC_ rl 

84, 



dC 



dC 

dip 



-d„ 



dC 



/ -i(Ttp a (x)^(y)) -i(Tijj a (x)iJjJ(y)) \ _ I S F (x,y) a/3 -iF(x,y) afS 
\-i( T il{x)ip{y)) -K T $l( x )$p(v)) ) \-iF(x,y) a p -S F (y,x)p a 

This is a 4x4 matrix. T means the time-ordered product, and (■■•) means the expecta- 
tion value. Sf is the Feynman propagator for quasiparticle, while — iF and — iF are the 
anomalous propagators. Next, we obtain the equations of motion for the propagators (13). 
We employ the Gor'kov factorization in (10) and (11), taking account of only the super- 
conducting pair-correlation by introducing the mean-field approximation. Then we obtain 
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the relativistically generalized (2+l)-dimensional Gor'kov equation written down as a 4x4 
matrix equation: 

i 7 ^ + 7 V A(x) \ I S F (x,y) -iF(x,y) 

A(x) t(rfd, - ( 7 °) T /i I \ -iF(x,y) S F (y,x) T 



(14) 



5^(x-y) 

5 ( V(x-y) 

A(x) and A(x) are 2x2 matrix mean fields, so called order parameters. The definitions are 

A(x ,x) a p = GF(x£,x;x ,x) a/ 3 = G(ip a (x^,x)ip^(xo,x)), (15) 
A(x , x) af3 = GF(x^, x; x , x) af 3 = G(4>l(x£, x)i>,3(xo, x)). (16) 

This gives the self-consistency condition. In general, the mean field clearly violates the 
Lorentz symmetry, as well as the gauge symmetry. In other words, the mean field involves 
quantities other than the scalar. 

We will also obtain the Fourier transform of the Gor'kov equation: 



% A \ I S F {k) -iF(k) 
A f J \ -iF(k) -S F (-k) T 



)-: 






: 



(17) 



Here k = (k + /i, k) and k = (k — fi, k). fk T means the transpose of ^. The self-consistency 
condition now becomes 

A = G J$r F ^ h - G Sik m (18) 

Here the mean field has only the internal degrees of freedom. 

The finite-temperature theory of the Matsubara formalism can be obtained in the same 
way. We introduce imaginary time r = it. The temperature Green's function is defined as 

g(x, y ) = -{T T 4?(x)*(y)) 

-{T T ^ a {x)i) P {y)) -(T T ^ a (x)^(y)) \ = / S(x,y) af3 -F(x,y) a p\ 
-(T T ^(x)4>(y)p) -(T T ^(x)^(y)) ) \-f(x,y) af3 -S{y,x)p a ) ' 

where (• • • ) means the statistical average. From the equation of motion of the temperature 
Green's function, the Gor'kov equation becomes 

-7°(|r -A*) +il k d k A{x) \ ( S(x,y) af3 -F(x,y) a/3 

A(x) "(7°) T (^ + fj) + ^{l k ) T ^ k ) { -?(x, yU S(V, x)fia 



(20) 



5 (3) (x-y) 

5^(x-y) 

Here the definition of the mean fields are the simple extension of those for the zero temper- 
ature: 

A(r, x) a/3 = GF(t + , x; t, x) = G(^(r+, x)$(t, x)), (21) 
A(r, x) Q/3 ee Of (r+, x; r, x) = G(^(r+, x)^(r, *)). (22) 

Fourier transform is also obtained as follows: 

7°(icj n + /i) - 7 • k A \ / <S(u; n , k) -F(uj n ,k) 

A (7°) T (^n - fi) - (7) T • * j V fc) -5(-w B , -k) T 



(23) 



1 
1 

Here (3 = (/c^T)" 1 ( k B ; the Boltzmann constant ), uo n = (2n + l)rr//3 is a fermion discrete 
frequency. Solving Eq. (23), we will obtain the solutions of (17) in the same form, except 
that we need to substitute ko — > icu n . 

IV. GROUP-THEORETICAL CONSIDERATION OF THE MEAN-FIELDS 

Now we perform the group-theoretical consideration of the mean fields. Under the Lorentz 
t r ansfor mat ion : 

V>V) = S^(x)=expC^Y)^), (24) 
then the mean field is transformed as 

<,\x')^ T <~' 

= SA(x)S 
i 
2 



A'(x') = (ij'(x')ij' T (x')) = (Sij(x)ij T (x)S T ) 



(1 + V 7 *)A(xXl + 



= A(x)+ l -u; u [r,A(x)C- 1 }C. (25) 
Thus we can decompose the mean fields as follows: 

A = (A 5 + A^)C, A = -C _1 (A 5 * + A^*7 M ), (26) 
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where S indicates scalar, while V indicates vector. Under the discrete transformations: 



ip(x ) ^ -i-f 2 ij(-x ) 



ip{xo,x) ^ -%P(xq,x')(-i^ 1 ), 
$(x ) ^ -4>(-x )(-i-f 2 ), 



(27) 
(28) 
(29) 



( here x = (xi,x 2 ) and x' = (—x 1 ,x 2 ) ), where C, V and T denote the operations of charge 
conjugation, spatial inversion and time reversal, respectively. Therefore, the mean fields are 
transformed as 



(^) ^ C(^ T >C _1 = -7 2 (^ T )7 2 
(^(x ,a;)^ T (xo,x)) ^(^(^o, x')if> T (x , x'))^ 1 
{ij T (x ,x)^(x ,x)) -7 1 (^ T (a;o,x')^(a;o,^'))7 1 



T\ ^ 2 



(30) 
(31) 
(32) 
(33) 
(34) 
(35) 



Thus each type of the mean fields is transformed under the spatial inversion and time reversal 

as 



r 



-AV, 
-A 5 * 7 2 , 

A V7 2 ^ A ^7°7 2 , 



r 



V„,l„,2 ^ 



Ar 7 x 7 



r 
r 



A ^7°7 2 , 

Ar*7S 2 , 



A l f A V 2 2 

2 7 7 -> A 2 V 7 7 , 



Ay* 2 2 
2 7 7 



(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 
(43) 



Therefore, with 2-dimensional rotation and parity, the mean field is decomposed into three 
irreducible representations: A 5 , A^ and (A]^,AJf). As expected, A s violates the parity 
similar to the Dirac mass term m^jjip. 
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V. THE GAP EQUATION 



In this section, we derive gap equations, and solve them numerically. For this purpose, 
first we have to solve the Gor'kov equation. Similar to the case of the (3+l)-dimensional 
theory [53-55], it is difficult to solve Eq. (17) or Eq. (23) completely in analytical form 
because of its matrix structure. Therefore we have to solve the equations assuming the type 
of the mean field that might be realized. Then we obtaine three Gor'kov equations for each 
type of the mean fields. These equations can be solved in the same way as the case of the 
(3+l)-dimensinal theory [53-55]. We give the following results. First the case of the scalar 
A s : 

' S F {k) -iF{k) 
-iF(k) -S F {-k) T 



1 / (g-|A s | 2 )^ A s Cfft- \A S \ 2 )C 



Next the case of Oth-component of vector Aq 



v . 



1 



/ 



D(k) 



S F (k) -iF{k) 
-iF{k) -S F {-kf 

- \K\ 2 )i°h° Kim - \k\ 2 i°)c 



(44) 



D{k) \-A s *C-\U-\A s \ 2 ) -C-\U-\A S \ 2 )}C 

D(k) = {k-k)(k-k)-2\A s \ 2 (k-k) + \A s \ A . (45) 



(46) 



^ -Al*C-\^ - K| 2 7°) -C- 1 ^ ^ - |A^| 2 )7 ^7 C 

D(k) = {k-k){k-k)-2\A^\ 2 {k-k + 2k 2 ) + \A^\\ (47) 
The case of lst-component of vector A\: 

S F (k) -iF(k) 
-iF{k) -S F {-k) rj 

i / CHh 1 - |AlThW An^ + |A^|V)c 



(48) 



D{k) \-&v*c-\fa^ + \&X\W) -c-'C^h 1 - |a^| 2 ) 7 1 ^7 1 c 

D(k) = {k-k){k-k)-2\AX\ 2 {k-k + 2kl) + \A\\\ (49) 
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In all cases, D(k) is second order in k%, and we can easily factorized it as D(k) = (k — 
E + )(k + E + )(k — E-)(k + EJ). Here, E + corresponds to the energy of the quasiparticles 
coming from the upper band ( positive energy states ), while E_ corresponds to the energy 
of the quasiparticles coming from the lower band ( negative energy states ). Because (44), 
(46) and (48) are 4x4 matrices, there is no degeneracy in these branches of the dispersion 
relations. This case relates to the fact that we treat a kind of spinless model. 

Now we construct the gap equations by using the Green's functions we have obtained. 
We use the finite-temperature Matsubara formalism. From the self-consistency conditions: 

VuvAO. (50) 



(3 J (2tt) 2 * 

we obtain the gap equation for the case of the scalar A 5 : 

G f d 2 k / 1 , /? 1 , /? \ 

1 = - / t — 777 ( -= =- tanh -E + + — — tanh -E_ J , 51 
2 J (2ir) 2 \2E + 2 + 2E_ 2 /' v ; 

E ± = v/(|fc| T/ i) 2 + |A^| 2 . (52) 

The second term in the integrand is the contribution coming from negative energy states 
and/or the lower band. In the context of this paper, the relativistic effect is the two-band 
effect. For the case of the Oth-component of vector A^, we get 

= G r d 2 k 1 {_ h (3 E +tanh d E _\ (53) 

2 J (2vr) 2 2v ^r^ypV 2 IF y J 

E ± = \k\ T J7+\Kf, (54) 
and the case of the lst-component of vector A] 7 , we get 

G f d 2 k f Jk 2 ^ + \AY\ 2 kl-kj ! s 



J fc 2 /i 2 + \AY\ 2 kj + k\ x q . 

+ ( V 1 =-)W~ tanh ^ E ~) ' (55) 



E ± = ^Jk 2 + fi 2 + \AY\ 2 =F ^/feV + lA^I 2 ^. 



(56) 



Integrals of Eqs. (51), (53) and (55) give positive quantities. Therefore, from these equations 
given above, we find there are possibilities to obtain the nontrivial solutions for the cases 
of the scalar A 5 and Oth-component of vector Ajf. To examine whether these equations 
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have nontrivial solutions or not, we have to check them more in detail. We will treat them 
numerically to study the characters of the solutions of these equations. We cannot obtain 
nontrivial solutions for the case of the spatial components of vector A\ and AJf, because 
the gap equations (55) become the form "— 1 = positive quantity". 
For the integration of our gap equations, we take 

r\ 2 h 1 r A r 27r 

! kdk de--- . (57) 



(2tt) 2 4tt 2 J Jo 

Therefore, our gap equations include four parameters: Coupling constant G, chemical po- 
tential /x, momentum cutoff A and temperature T. Because our theory treats a massless 
case, there is no unit of energy in our theory. To see the effect of the lower band, we also 
treat the gap equation of the "no sea" ( neglect the contribution of the lower band (Dirac 
sea) in the gap equation of the scalar ) case: 

G f A d 2 k 1 , 8 

E + = y/(\k\-n) 2 + (A™™ 1 ! 2 . (59) 

There is a simple scaling relation in our gap equations. When we transform the gap 
equation in next relations: 

\k\/n=\k>\, A//i = A', j3iM = j3', |A|//i=|A'|, (60) 

Eq.(51) is transformed as 

1 = —fi f k'dk'i ; 1 tanh-./nfc 7 ! 

+ ; --, 1 tanh^J(|fc'| + l)2+|A'| 2 V (61) 

y/(\k'\ + 1) 2 + |A'| 2 2V VI i ; 

Therefore, when the ratio A/// is fixed and when we take G' = G/i = const., we treat the 
same equation with (61). 

Usually, a gap is much smaller than the Fermi energy: |A| <^ fi. For example, in the 
case of solid, the ratio |A(T = 0)|//x is 10~ 4 — 10~ 2 . We have to choose model parameters 
G and A to satisfy the condition. In principle, we can choose A up to the energy scale 
where the conical dispersion is a good approximation. Both the Fermi momentum fop and 
the upper bound of A should be in the order of the inverse of the lattice constant cr 1 . Here, 
we choose the momentum cutoff as A = 2kp = 2/x, and fix it throughout this study. After 
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k /|2 



choosing the values of parameters, the integration in the gap equation is performed, and 
search the self-consistency condition under the variation with respect to the amplitude of 
the gap. We have performed the integration in our gap equations by using the numerical 
package Mathematica version 4.1. 

Figure 3 shows the gap at T = as a function of the coupling constant G. Here we take 
the momentum cutoff A as A/2// = 1 and we set /i — 1. Similar to the usual nonrelativistic 
BCS theory, both the scalar and "no sea" depend exponentially on G. The gap of the scalar 
is always larger than the "no sea" case. The ratio \A S (T = 0)|/|A"°- sea (T = 0)| is almost 
always 1.57. 

In Fig. 4, we show the temperature T dependence of the solutions in the scalar and "no 
sea" cases. Here we set G = 2/i -1 , /t — 1 and A/2// = 1. In the cases of graphite or MgB 2 , 
the relativistic model can be applied to the region of the energy width almost 1 — 2eV. If 
we take /x = leV ( A = 2eV ), the critical temperature T c becomes 38K (0.0033eV) for 
the scalar and 24K (0.0021eV) for the "no sea". The ratio T° calar /T™~ sea is 1.57. Both 
the scalar and "no sea" cases fulfill the BCS universal constant |A(T = 0)|/T C = 1.76 [57]. 
Therefore, both cases obey the BCS-like temperature dependence. 

From our numerical results, we conclude that, the scalar ( the two-band case ) gives 
always larger solution than the "no sea" ( the one-band case ). The contribution coming 
from the lower band ( the negative energy states ) enhances the superconducting gap. It 
should be noted that, the values of the amplitude of the gap and the critical temperature 
do not predict the superconductivity in real substances directly, because it is impossible to 
remove the arbitrariness of the choice of the values of the model parameters completely. 
However, when we consider the two-band system, we have to keep in mind the fact that the 
lower band gives a sizable effect in the superconductivity. 

About the Oth-component of vector case, we could not find a reasonable solution in our 
numerical calculation. The reason of this result should be understood in the following way. 
We rewrite Eq. (53) at T c as follows: 



The integral of the right hand side gives the conserved charge of the system. Therefore, the 
equation for determination of T c has no temperature dependence: We cannot determine T c 




(62) 
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by Eq.(62). In the (3+l)-dimensional massive theory, the Oth-component of vector pairing 
gives a meaningful state, and we discussed various aspects of the relation between the scalar, 
"no sea", "nonrelativistic" and Oth-component of vector [53-55]. This fact contrast with the 
results of the (2+l)-dimensional theory. 

VI. THE KOHN-LUTTINGER EFFECT 

Now we consider the pairing problem with the repulsive interaction G < in the La- 
grangian (2). For this problem, we should examine the appearance of a pole in the 4-point 
function ( the 2-particle Green's function, see Fig. 5(a) ) by using the Bethe-Salpeter ( BS ) 
formalism [44,58], like the work of Kohn and Luttinger [45]. We start with introducing the 
fermion-fermion BS equation [59] in the finite-temperature Matsubara formalism: 



We treat (63) as a function of temperature T, and the critical temperature T c is deter- 
mined when a self-consistent solution of x{p) appear in Eq.(63). The normal-state fermion 
propagators are given as follows: 

S(u n ,p) = -I, S(-u n , - P ) T = -C- X \c. (64) 
V V 

Here p = (iu n + /i, p) and p = (iuo n — //, p). T(p,p') is the irreducible vertex part. The BS 

amplitude x(p) * s decomposed as follows: 

4 

X(P) = (X S (P) + X V ,(phnC = (J2 XaT A )C. (65) 

A=l 

Here, S and V denote scalar and vector, respectively. x(p) nas to fulfill the Pauli principle: 




(63) 



Xpa(-P)- 



(66) 



Then, each component should obey the following relations: 



X S {p)C aP = -x s {-p)c?, 



x;(p)(rcu 



x;(-P)(rcw 



(67) 



Therefore we get 



x s (-p) = x s (p), 



xl(-p) = -xl(p)- 



(68) 
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In Eq. (63), replace x(p) to the expanded form ( Eq. (65) ), and take trace of both sides, 
we obtain the following form: 

XA{ P )\tr{Y A Y A ) = ^f( P , P ')^tr(r^K,p')XA(p / )r A 5(-^,- P ') T )(69) 

Similar to the treatment of the Gor'kov equation in Sec.V, we completely neglect couplings 
between different types of pairing functions. Then, for specific A, we obtain the following 
equations: 

x (25^0 tanh * w ~ " ] + tanh ^ m + ^ ' (70) 

X V o{p) = \ J ^f(p,pOXo y (p')^(-tanh^(|p'|-^)+tanh|(|p'| +/ i)), (71) 



Mp'\ 



In this paper, we estimate t(p,p') by using a random phase approximation ( RPA ) ( as 
illustrated in Fig. 5(b) ). The polarization is given as 

1 f 

n(wi,9) = ^2^J j^trS(ui + u n ,p + q)S(u n ,p). (73) 

Here, lo\ is a boson discrete frequency. After summing up cu n , we negrect u>i dependence of 
TI(uji, q) to obtain the static polarization: 

f d 2 p 1 



11(0, q) 



x 



(2vr) 2 6 p+q -6 p 

p^L( I + I A 

e p \ e P( € p-ri + 1 e P(tp+ri + i J 

{P + Q) Q ( 1 + 1 A \ 

e p+q \ e ^p+q-n) + 1 e fcp+q+ti + l )y 



(74) 



where, ep = \p\ is the relativistic dispersion for massless particle. The integration in (74) has 
a ultraviolet divergence, and we have to introduce a momentum cutoff A. To obtain 11(0, q), 
we perform the integration numerically in Eq. (74). Figure 6 shows the results of numerical 
integration for 11(0, \q\) with several T. Here we set A = 2p F = 2/j, with ji — 1. We find 
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a peak near \q\ = 2p F . This behavior is a reflection of the sharpness of the Fermi surface. 
The peak decreases when T increases. To treat the Cooper problem, we concentrate on the 
behavior of 11(0, \ q\) at the Fermi surface. To see this, we define q = p' — p and substitute 
\q\ to \p' — p\ in n(0, |g|). Then we set \p\ = \p'\ = p F , we get U(0,p F ,PF, cos ^p.p')- 
Here Qp.p' = @p ~ @p'- We calculate U(0,pf,Pf, cosOpp) numerically. Figure 7 shows the 
angular dependence of U(0,pf,Pf, cosdp.p') with several T. The dent at Qp,p> = vr grows 
when T incleases. At enough low temperatures, the angular dependence is almost cos 9, and 
then we can use the following expression to examine the pairing properties of the system: 

U{0,p F ,p F ,cos9pp') T=0 ~ ^-(A - p F ) + a(l- cos 9pp>) 

= — (A — p F ) + a — a(cos9p cos9pi + sin 9p sin 9 pi). (75) 

From Fig. 7, we find a ~ 0.0284. Then we obtain the expression for T(p,p') at low 
temperatures within the RPA: 

G 

T{p,p,) ~* i-Gn(o, PF ,p F ,cos0 p . pl y (76) 

To obtain the BS equation for specific symmetry of pairing, we employ the angular decom- 
posion to x(p) an d ^(p, p'). The two-dimensional angular momentum eigenfunction is given 

as 

Vi(0) = -^f l \ (77) 



where, / is an angular momentum quantum number in a two-dimensional system. We only 
consider the time-reversal invariant pairings because of simplicity. Therefore, we take a 
linear combination: 

' 1 '° + e~ m ) = cos 19. (78) 



Each component of the amplitude is decomposed by this function: 

X S (P) = E X S (\p\)i^co S l9 p -> £ xf^cosWp, (79) 

l-.even b.even 

Xo(p) = J2^(\p\)i^^sl9p^Y,^o)i^^ S l9 p , (80) 

l-.odd hodd 

\X V {P)\ = El^(|Pl)l^ C0S ^P^E^ y |^ C0S ^P' (j = Z± 1). (81) 

v 71 " v 71 " 

y.oda j -.odd 
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Here, we negrect \p\ dependence of the amplitudes for weak coupling approximation: 
IxOpDI* ~ \x\i- Because of (68), we choose even / for scalar in Eq.(79), while we choose odd 
I (j) for vector pairings in Eqs.(80) and (81). The irreducible vertex part is also decomposed 
as follows: 

1 1 



G- 1 - (A - p F )/2ir - a + a cos Opp G' 1 - (A - p F ) /2tt - a 

a 



(G- 1 


-(A 


-p F )/2ix 


-a) 2 






a 2 




' (G- 1 


"(A 


-Pf)/2tt 


-a) 3 






a 3 




(G- 1 


"(A 


-Pf)/2tt 


-a) 4 



cos9 pp/ 
cos 2 9p p 
cos 3 Op p 

+ •••• (82) 

Therefore, when we decide the quantum number / ( or j ) of a BS amplitude, the components 
coupled with the effective interaction T(p, p') are determined a priori : From the form of the 
expansion given in Eq. (82), we find the fact that, 1st, 3rd,... terms in Eq. (82) couple only 
with components of even /, while 2nd, 4th,... terms in Eq. (82) couple with components of 
odd It is an important fact that, the sign in each term in the expansion for T alternates 
between plus and minus. 

Based on the preparation given above, we obtain the BS equations for specific / ( or j ) 
of several types of pairings. For the case of the scalar: 

1 = 



x 



(2(E?ho tanl ^ (|p ' 1 " + wr^ tanh f M + "')- (83) 



For the Oth-component of vector: 

r-A 



1 8tt 2 



h*f P> , ^(- tanh f(IP'|-A i )+ tanh f(IP'l ( 84 ) 



Here / is determined as 



( 2 *n COBZfljCOBZfla 

d °i / ^27^ 7T ~ Wo _ — — — ~ Q , (85) 



^0 



G 1 — (A — pf)/2h — a + a cos 6*12 ' 



where 0\ = Op, O2 = 0p> and 0\i = 0\ — 2 . In principle, we have to choose A in our 
BS equations (83) and (84) as the same value for the polarization in Eq.(74). By using 
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the expansion (82), when G < 0, we find the following fact: In the angular integration of 
/, there are only negative contribution for I ( I < ) of the scalar case, while there are 
only positive contribution for I ( I > ) of the vector case. Therefore, from Eq.(83), we 
recognize that there is no solution in the scalar pairing. From the same reason given in the 
discussion for Eq. (62), we cannot find solution for the Oth-component of vector case. For 
the spatial-component of vector, we get 

Wo P ' dP 'l * 

cos j9i cos j9 2 



G- 1 - (A - p F ) /2tt - a + a cos 9 



12 



p'\ 2 cos 2 9 2 — fi\p'\ , P ,, \p'\ cos 2 9 2 + fJ,\p'\ P 



tanh i (|p1 - " } - 1 TpV;; tmh i (ip,|+ ^ 

In Eq. (86), the angular integration gives the negative quantity, but the integrand of mo- 
mentum integration is also always negative with respect to the variation of T. Therefore, 
it is impossible to find T c , and we have no solution in the spatial-component of vector. 
The summary of the results of the solutions of various types of our gap equations and BS 
equations are presented in table I. Only the scalar pairing in the attractive interaction can 
have nontrivial solution. In our treatment ( BS-RPA ), we could not find a KL effect in 
our model. Baranov et al. [60] obtained a conclusion that no pairing instability arises in 
a two-space-dimensional nonrelativistic model with second order perturbation. Chubukov 
performed a calculation beyond second order: He included a vertex correction, and showed 
that the KL effect arises in the two-space-dimensional Fermi liquid [47] . Our negative result 
of the KL effect is based on both the single conical dispersion model and the angular depen- 
dence of the irredicible vertex part estimated by RPA. Therefore, our result does not deny 
the presence of the KL effect in (2+l)-dimensional systems completely. In fact, Gonzalez et 
al. showed the presence of the KL effect in graphite, by taking into account the interaction 
between two inequivalent Fermi points carefully [32]. On the other hand, as discussed in the 
introduction, we have set the purpose of this paper just to examine the lower-band effect in 
superconductivity. For the purpose, we use a single relativistic-dispersion model with [i ^ 0. 
The situation of our model is different from that of the work of Gonzalez et al. 
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VII. CONCLUDING REMARKS 



In this paper, we have investigated the relativistic model of two-band superconductiv- 
ity in (2+l)-dimension. After we have introduced the model Lagrangian, we have derived 
the (2+l)-dimensional relativistic Gor'kov equation for the pairing problem under the four- 
fermion contact attractive interaction. We have performed the group-theoretical considera- 
tion of the gap functions. The characteristic aspects of the gap functions were revealed. By 
using the Gor'kov formalism, we have derived the gap equations for several types of pairings, 
and have solved them numerically. We have understood quantitatively the effect of the lower 
band in the two-band superconductivity: We have found the lower band enhances the su- 
perconductivity. We also have examined the pairing problem under the repulsive interaction 
( the KL effect ), by using the BS formalism. 

Now, we discuss some remaining problems and/or further possible investigations in our 
theory. It is interesting to perform the calculations of the response function or polariza- 
tion function under the presence of electromagnetic field, because of the reason that those 
functions can give the Chern-Simons ( CS ) term or not. Goryo and Ishikawa discussed the 
induction of the CS term in (2+l)-dimensional nonrelativistic theory, with parity and time- 
reversal violating superconductors [61]. Hosotani showed the (2+1) -dimensional QED with 
the CS term dynamically generate a magenetic field [62]. Miransky et al. studied the fact 
that an external magnetic field enhances a fermion dynamical mass, and this phenomenon 
is universal in any models of (2+1) and (3+1) dimensional field theories [63]. We suppose 
these studies should have an intrinsic relation. We have a plan to investigate these physics 
in our model in near future. 
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FIG. 1: (a) The structure of the honeycomb lattice. A and B denote the two different sublattice 
sites, (b) A schematic figure of the band structure of the honeycomb lattice. By symmetry, the 
hexagonal two-dimensional Brillouin zone of it has two degeneracy points. In the case of graphite, 
the Fermi levels locate on these points. In this paper, we use a Dirac fermion model with /t / 0. 
Our theory can be applied to electron-doped graphite, etc. 

FIG. 2: A schematic figure of the relativistic dispersion in a two-dimensional system, ep denotes 
the Fermi energy of the system. 

FIG. 3: The G dependence of the pairing gap at T = 0. We set A = 2k f = 2/t and // = 1. 

FIG. 4: The temperature dependence of the pairing gap. We set G = 2 / u~ 1 , /t = 1 and A/2// = 1. 

FIG. 5: (a) The diagrammatic representation for the two-particle Green's function G^ . f denotes 
the irreducible vertex part, (b) The diagram for the irreducible vertex part within the RPA. The 
solid points represent the bare vertices G. 

FIG. 6: The (/-dependence of the static polarization 11(0, q) under various temperature. Here we 
set /t = 1 and A/2// = 1. 

FIG. 7: The ^-dependence of the static polarization 11(0, kp,cos6) under various temperature. 
Here we set \x = 1 and A/2// = 1. 
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